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1 INTRODUCTION 



The very first stars in universe were borned in the lack of 
heavy elements. It could leads to born of the super mas- 
Tl ' sive stars (SMS) (iBromm et al.ll20ul; iBromm fc Loebll200j ; 
> I ISantoro fc Shull I l200rj ; iDevecchi fc Volonteril |2009h . Such 



stars can be as massive as 1000M Q < M < 5OOOOM . 

The stars with masses about IOOOA/q, which will be 
referred to as Very Massive Stars (VMSs) are expected 
to collapse into black holes with very little mass loss 
l|Frver et alj|200ll ). They would leave behind massive black 
holes (MBHs) , which could play the role of seeds for the su- 
per massive black holes (SMBHs) of Active Galactic Nuclei 
(AGNs ). The collapse of VMS were d iscussed recently in the 
paper l|Komissarov fc Barkov lliulOh . 

Even more massive 3 x 1O 4 M < M < 10 6 M© stars, 
which will be referred to as Super Massive Stars (SMSs) 
could be formed in more massive dark matte r haloes with to- 
tal mass M ~ 10 s Mp, collap sed at z ~ 10 l)Bromm fc Loebl 
120031 ; iBegelman et all 12009 ). This stars do not reach the 



ABSTRACT 

The very first stars in the Universe can be very massive, up to 10 6 Mq. They would 
leave behind massive black holes that could act as seeds for growing super massive 
black holes of active galactic nuclei. Given the anticipated fast rotation such stars 
would end their live as super massive collapsars and drive powerful magnetically- 
dominated jets. In this paper we investigate the possibility of observing the bursts of 
high-energy emission similar to the Long Gamma Ray Bursts associated with normal 
collapsars. We show that during the collapse of supercollapsars, the Blandford-Znajek 
mechanism can produce jets as powerful as fewxl0 51 erg/s and release up to 10 56 erg 
of the black hole rotational energy. Due to the higher intrinsic time scale and higher 
redshift the initial bright phase of the burst can last for about 10 5 seconds whereas 
the central engine would remain active for about 10 days. Due to the high redshift the 
burst spectrum is expected to be soft, with the spectral energy distribution peaking 
at around 60keV. The peak total flux density is relatively low, fewxl0 _7 ergcm _2 s _1 , 
but not prohibitive. The such events should be rear 0.03 year -1 , the observations needs 
long term program and could be done in future. 

ISuwaet al.l[2007h . The other is that the accretion disks 
of supercollapsars are far too large and cool for the neu- 
trino annihilation mechanism. This has already been seen 
in the numeric al simulations of s upercollapsar with mass 
M = 300M© dFrver et alj l200ll). U t ilising the study of 
hyper-accreting disks by Beloborodovl (I 20081 ) we find that 
at best the rate of heating due to this mechanism is 
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6 x 10 43 M 9/4 M-6 /2 erg s ~\ 
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of SMS can provide an alternative way of producing SMBHs. 

There are two crucial differences between a normal col- 
lapsar and a supercollapsar. One is that instead of a proto- 
neutron star of solar mass the supercollapsars develop proto- 
black holes of tens of solar masses, within which the neutri- 
nos from electron capture are trapped (|Frver et al.l l200ll : 



where M is the accretion rate and Mh is the black hole 
mass. (Here and in other numerical estimates below we use 
the following notation: is the mass accretion rate mea- 
sured in the units of lO fe M0S _1 and Mk is the mass mea- 
sured in the unit s of lO fe M0.) Such low values have lead 
iFrver et al.l ()200ll ) to conclude that the magnetic mechanism 
is the only candidate for producing GRB jets from supercol- 
lapsars. In the following we analyse one particular version of 
the mechanism where the jets are powered by the rotational 
energy of the black hole v i a the Blandford-Znajek proc ess 
iBlandford fc Znaiek|[l977l : iBarkov fc Komissarov IboOSl ). 

This paper is an ext ension our previous work 
l|Komissarov fc Barkov I l20ld) here we investigate proper- 
ties of collapse of VMS. In this paper we estimate temporal 
structure of the SMS collapse and predict the observational 
evidence of this process. 
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2 PHYSICAL MODEL 

Lets look for this problem. SMS should be fully convective 
and have rigid rotation. In our simple calculations we neglect 
aspherical shape of fast rotation SMS. 

We assume the radius SMS before collapse as 640 
GM/c 2 dZeldovich fc Novikovl Il97ll : IShapiro fc Teukolskvl 
Il983l ; IShibata fc Shapirol |2002h . We can derive the undi- 
mensional critical parameter of rotation for SMS with 3 x 
1O 4 M < M < 1O 6 M as J a < (M/lO 6 Af ) 1/4 . Numerical 
simulations of SMS collapse show forma tion of massive BH 
and accretion disk . Follow the w orks |Shibata fc Shapirol 
2002; iBarkov fc KomissarovlboiOh we know spin parameter 
of BH and mass of accretion disk. 

In sake of simplicity we use model of the advecti on dom- 
inated accretion flow (ADAF) (|Naravan fc Yill 19941) . 

Assuming gas pressure to dominate, we write the pres- 
sure P = pc 2 where p(R) is the height-averaged density and 
c s (R) is the isothermal sound speed, we also put 7 = 4/3. 
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R ' 
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2 2 2 



(2) 



here v r is radia l velo city, a is a constant 
l|Shakura fc Sunvaev I I 19731 ). Vertical scale height is 
H ~ Rc s /vk- It is easy to see that disk accretion time can 
be expressed as 
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14 / R d 



(3) 



3v r 9a V GM a 

Lets assume accretion disk mass as Md = DM S , the accre- 
tion rate is determined by defuse time from maximum disk 
mass radius 



t d 



14 V R 



(4) 



Using continuity equation M = —4nRHv r p and exp. Q, 
get 



(5) 



_ 3V14 DM S 
P ~ 16vr R3/2R^ 2 ' 

combine ([2]) and ([5]) we get 



3V14 DGMj 

56n 7^5/2^3/2 



(6) 



magnetic field pressure shou ld be the same order as gas one 
B 2 /8tt w P, follow the work lTout fc Pringle I (jl996l ) we can 
estimate large scale poloidal magnetic field B v m B C H/R, 
so 

3 l/2 2 3/4 ^Q Ms 
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73/4 R5 /4 R 



3/4- 



After last marginally stable orbit the large scale mag- 
netic field could be accreted to BH with ga s from accretion 
disk, the similar approach are described in iRevnolds et al.l 
l|2006h . 

Rms is the ra dius of the marginally stable orbit 
(|Bardeen et aUll972ft . 



Rms = r g {3 + Z 2 ^ [(3 - Zi)(3 + Zx + 2Z 2 )] 1/2 }, 

Zy = 1 + (1 - a 2 /A^" 2 ) 1/3 [(l + a/M) 1/3 + (1 - a/M) 1/3 ] 

Z 2 = (3a 2 /M 2 + Z 2 ) 1/2 . 



(8) 



To estimate magnetic flux on the horizon of BH we 
should substitute JSJ to ([7]) so 
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3 1 / 2 2 3 / 4 R 3 JfVDGM s 

7 3/4 ^3/4 



(9) 



t he energy realise of Blandford-Znajek mechanism 
f l|Blandford fc Znaiek|[l977h ') could be expressed as 



(10) 



here fl = flh/2 the angular velocity of field lines, where 
O fe = a/2(l + VI - a 2 )c 3 /GM B H is an gular velocity of BH 
( Komissar ovl |200"3 ) . a = Jbh / M% H is dimensionless speen 
parameter of BH. Let us express the current radius and disk 
radius in r g units as r ms — Rm S /r g and Td = Rd/r g , substi- 
tuting |9j to JTDJ we get 



Ebz 



2 3/2 



D 



3136 (1 + Vl - ^) 2 G r 3/ 2 



(11) 



The energy release depend on spin parameter of BH a, mass 
fraction of the disk and radius of the disk. It is easy to see 



Hard X-ray Bursts from Collapse of the Super Massive Stars 3 



that all dependence from 0.5 < a < 1 is very weak (0.74-1.4) 
so for our continue estimations we can assume it as 1. The 
expression will be shown as 



as follow r) = M/Mi n it « (r in /r ou t) p , here r in « r„ 



5r g , 



Ebz 



c D 

0.0012 



(12) 



3 RESULTS 

We use the method which is described in the paper 
iBarkov fc Komissarovl ((2010) to calculate accretion disc 
fraction and accretion time scales. The accretion disk which 
formed for spin parameters of the star 0.5 < J s < 0.97 have 
longer life time (see fig. [T] left panel) then free-fall time 



// 



2gj 
9GM S 



(13) 



Parameter D depends on star angular momentum. We in- 
terpolate our numerical model for polytropes (7 = 4/3) star 
with good accuracy in the range 0.5 < J s < 0.97 by simple 
formula 



,D = -3X + 7J S - 4.9 



(14) 



. For J s < 0.5 the accretion disk mass is negligible and 
£><1. 

The accretion disk has maximum of matter distribution 
near radius Rd ~ 25 ~ 50r 9 and it is a wick function of 
J a (see fig[T] right panel). To simplify our estimations we 
put Rd ~ 40r s But accretion could be 3 4- 10 times longer 
(see fig[T]left panel) or you can find scaling for self similar 
solution in the work of iBeloborodovl (| 20081 ). 

The disk accretion time we can estimate as 

t d ~ lM Si6 a_i _1 day (15) 

Using the mass accretion rate of fig. {TJ we can check 
if the neutrino cooling needs to be included in the model. 
Under the conditions of the supercollapsar's disk its cooling 
is dominated by pairs. Using the well kn own equation for this 
cooling rate (e.g. lYakovlev et al. II2001T I we can compare the 
cooling time with the accretion time at a given radius. The 
result is 



t d 



IQ-^orJ/^R/Rg) 



-13/8^5/4 M - e 



3/2 



(16) 



Thus, except for the very inner part of the disk, the neutrino 
cooling is indeed inefficient. 

The complicated magnetic field topology can leads to 
strong degradation of energy release, the ene rgy flux could 
be 10 times lower (tBarkov fc Bau shcv 2009). Over factor 
which can decrease magnetic flux on horizon is runaway o f 
magnetic field if it became to strong l|lgumenshchev Il2008h - 
We will introduce the magnetization parameter in our for- 
mula P which could be order 0.1. 

Here, we assumed that the whole of the disk is accreted 
by BH, following the original ADAF model. However, it has 
been argued that this model has to be modified via including 
the disk wind (ADIO S (Advection Dominated Inflo w Out- 
flow Solution) model, Blan dford fc Begelman|[l999l ). which 
implies a mass loss from the disk and a smaller accretion 
rate compared to FigfTJ While the arguments for disk wind 
are very convincing, the actual value of mass loss is not well 
constrained. We can estimate the effectiveness of accretion 



p « 0.75. In the end we have 77 ~ 0.1. 

Using (112[) now we can estimate The maximum posssible 
energy realise rate is 

Ebz ~ 1.7 x 10 51 7/_i/3_iD_ierff s~\ (17) 

here £>_i = D/0.1, 77-1 = 7//0.1 and /3_i = /3/0.1. The 
remarkable property of disk dynamo supported BZ is unde- 
pendable from mass of the progenitor. It is mostly depend 
only from SMS angular momentum (|14l) . Higher angular mo- 
mentum - longer and brighter event (see fig[T]left panel). 

The total power of such event is enormous 
Etot ~ td x Ebz ~ 10 56 ergs. This energy is enough to 
sweep away the protogalaxy with baryonic mass ~ 10 s Mq 

Taking in to account jet collimation, we can expect very 
bright sources, which could be detected from very large cos- 
mological distance. 



4 DISCUSSION 

SMS can be borne only in very early stages of star formation 
[z ~ 10), the observational time should be 



tdo = (1 + z)td « 10 days 



(18) 



Due to large z the spectra should be shifted to Hard 
X-Ray. If our analogy with short-hard GRB is right then we 
expect spectral peak in aria 60keV 

We do not know for sure what is the central engine for 
short-hard gamma ray burst (GRB), but it is most probable 
is merging of two compact stars, accretion disk formation 
and collimation of the jet produced by strong wind from the 
disk. Here we may expect the similar behaviour. The same 
collimation properties leads us to opening angle 0.1 rad. Put 
the factor dilution A ~ 10 -3 

Let us estimate brightness of such burst. For [z ~ 10) 
and Cosmological parameters Op = Oa + f2 m = 1 and J1a = 
0.72, Q A = 0.28 jKomatsu et aljr2 009) lumi nosity distance 
will be d L » 3.2 x 10 29 cm |Mukhanov Il2005l 1. so flux on the 
earth should be 



F = 



Ebz 
\-xd\A 



1.3 x 10" Vie c 



i/3_iD_ierg cm 2 s 1 



(19) 



here e c ,-i = e c /0.1 is conversion coefficient of Blandford- 
Znaek flux to radiation. SMS were for med first in gas 
clouds with mass ~ 10 s M Q . Follow the wor k lBegel man et al.l 
(2006) we can estimate the density of SMS at z = 10 as 
n mh ~ O.OlMpc" 3 . 

We can estimate the rate of su ch events in analogy to 
work iKomissarov fc Barkov I (|2010l ). Let us assume, for the 
sake of simplicity, that all supercollapsars go off simultane- 
ously at cosmological time t e corresponding to z = 20 (a 
moderate spread around this redshift will not significantly 
change the result). In flat Universe the observed time sep- 
aration between events occurring simultaneously at ro and 
r\) + dro, where ro is the commoving radial coordinate, is 
dt — cdro- The corresponding physical volume within one 
steradian of the BAT's field of view, is 

dV = a 3 (t e )rldr , 

where a(t e ) = (1 + z)^ 1 is the scaling factor of the Universe 
at t = t e (in the calculations we fix the scaling factor via 
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the condition a(t ) = 1). r and t e are related via r = 
r L (l + z) -1 - Putting all this together we find the rate to be 

f c = A^,0.0SA-s^. (20) 

This events should rear and the long term observations 
program could reveal them. 



5 CONCLUSIONS 

In spite of the significant progress in the astrophysics of 
Gamma Ray Bursts, both observational and theoretical, it 
may still take quite a while before we fully understand both 
the physics of the bursts and the nature of their progenitors. 
At the moment there are several competing theories and 
too many unknowns. Similarly, we know very little about 
the star formation in the early Universe. For this reason, 
the analysis presented above is rather speculative and the 
numbers it yields are not very reliable. Further efforts are 
required to develop a proper theory of supercollapsars and 
to make firm conclusions on their observational impact. The 
collapse of SMS could be one of the most powerful event 
in the universe. This event can destroy the seed cluster 
and form single SMBH. The expected very long duration 
of bursts and their relatively low brightness imply that a 
dedicated search program using the image trigger will be 
required. Such search would be useful even in the case of 
non-detection as this would put important constraint on 
the models of early star formation, GRB progenitors, and 
SMBHs. 
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